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maximally entangled stabilizer states as MUB vectors. Consequently, these bipartite entangled stabilizer MUBs 
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I. INTRODUCTION. 

One of the most important and long-studied tools in quan- 
tum information theory is that of mutually unbiased bases 
(MUBs). Two orthonormal bases A = {\a}} and B = {\b)} 
in a Hilbert space of dimension d are said to be mutually un- 
biased when |(a|fe)| = 1/V" i-e. certainty of a measurement 
outcome in one basis implies complete uncertainty of a mea- 
surement outcome in another. This is the finite-dimensional 
analogue to the complementarity of position and momentum 
in continuous variable quantum mechanics. Typically, MUBs 
are most useful in Hilbert spaces, W c /, of prime power dimen- 
sion (d = p k ), for which complete sets of MUBs are known 
to exist and a number of construction methods are available. 
Ignoring the trace component (which is often known or unim- 
portant), decomposing a d x d Hermitean operator (e.g. a 
density matrix) requires d 2 — 1 parameters, which necessi- 
tates measuring d + 1 different observables (since each ob- 
servable yields d — 1 independent probabilities). Complete 
sets of MUBs are sets with d + 1 orthonormal bases, possess- 
ing the desirable properties of being both mutually unbiased 
with respect to one another and also being minimal in terms of 
the number of observables required (hence this is considered 
the optimal tomography set-up 01 lit)- 

Our work here concerns the construction of MUBs in 
Hilbert space of dimension d~p 2 that are deliberately incom- 
plete in that they contain only p 2 — 1 observables - insufficient 
for parameterizing all operators in %„2, but sufficient and min- 
imal for the description of Hermitean operators that are local 
maximally mixed (LMM) y|]. LMM operators, W, defined 
on a bipartite system H p 2 = CP ® C p are those for which 
Tri(W) = Tri(W) oc I. This class of operators is surprisingly 
broad. The Jamiolkowski isomorphism, for example, tells us 
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that any unital map £ acting on H p can be represented by 
an LMM operator, indicating that this result could potentially 
be useful for the characterization of noise processes, whilst 
reducing the number of measurements required (process to- 
mography using a similar construction is discussed in detail 
in 101). Other scenarios in which the non-local information 
is of paramount importance include investigation of bipartite 
entangled and non-local states, and the witnesses lUt] and Bell 
inequalities [6Q that identify them. As a final example, the mo- 
tivation for this work came in considering a convenient, min- 
imal basis with which to decompose so-called Clifford wit- 
nesses for detecting stabilizer vs. nonstabilizer operations 10]. 

The literature concerning MUBs, constructions and related 
structures is vast. This field of study seems to originate with 
Schwinger's construction for unitary operator bases (0 in 
1960, and subsequently Ivonovic's 1981 construction lUOll for 
complete MUBs in prime dimensions. Wootters and Fields 
provided a MUB construction for power-of-prime dimen- 
sions d = p k and showed its optimality for state reconstruction 
(tomography). A more recent (2002) construction that also 
works for power-of-prime dimensions is given by Bandy opad- 
hyay et al. Hill and this is framed explicitly in terms of Pauli 
operators and stabilizer states. Lawrence et al. Il2ll found a 
similar construction for multi-qubit systems in the same year. 
Since then a large number of related results have been pub- 
lished e.g. HH^T^] (also see a recent review article lfl7ll and 
references therein) and a number of interesting connections 
with combinatorics (e.g. mutually orthogonal Latin squares 
lfl8|p and finite geometry lfl9Tl2lll . A prominent example of 
the usefulness of MUBs is their optimality for state or pro- 
cess reconstruction (a recent experimental result tfj]] shows an 
improvement over standard techniques by using MUB state 
tomography). Quantum key distribution schemes II27[ 12811 typ- 
ically rely on MUBs for their security. Another important 
application of MUBs is their interpretation in terms of finite 
phase space, leading to a discrete Wigner function; for a par- 
ticular choice of MUB using stabilizer states, the resulting 
Wigner function can shed light on the computational power 
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of circuits in the so-called "Clifford computer" model 122V 
l2l. 



H. DEFINITIONS AND USEFUL RESULTS 
A. Relevant finite groups and their properties 



Inadvertently, we have rediscovered some results that were 
prev iously known in the context of quantum key distribution 
ll28ll . and in the context of unitary designs 1RI I29I1 (i.e., the 
Cliffords that we use to create some of our BES-MUBs are 
known to create a minimal unitary design). Recent work by 
Planat [30] is somewhat related to our current investigation, 
insofar as it utilizes graph theoretical concepts and stabilizer 
(Pauli operator) observables to examine the construction of 
MUBs. Kalev et al. lfl6ll investigated MUBs in bipartite sys- 
tems using sets of commuting Pauli operators, but their work 
is more focused on complete sets of MUBs for density opera- 
tors in V. 1. 



This work provides an alternative graph-theoretic method 
(as opposed to unitary designs or finite field constructions) of 
analyzing MUBs and similar structures in quantum informa- 
tion theory. It is hoped that a combination of the alternative 
methods outlined here, in addition to those of I28l430ll and 
others, will prove fruitful for further analyses. We show how 
to create an orthonormal basis of p 2 stabilizer states in ~H.pi, 
given a matrix F € SL(2,Z P ). Furthermore, we show that the 
quantity Tr(-F. _1 F;) indicates whether the bases corresponding 
to Fi and Fj are mutually unbiased. This leads naturally to a 
Cayley graph structure wherein graph vertices are given by the 
elements of SL(2,Z p ), and edges between vertices correspond 
to mutual unbiasedness of the corresponding bases. The BES 
MUBs that we seek are easily shown to be maximum cliques 
of the Cayley graphs, and for primes up to 11 we can parti- 
tion SL(2,Z P ) into p distinct (non-overlapping) BES-MUBs. 
For primes 13 and higher, it is an interesting open question 
whether such BES-MUBs exist, as a deterministic search for 
the maximum clique is infeasible. For the related question of 
minimal unitary designs it has been noted by Chau that sub- 
groups of SL(2,Z p ) of a particular size only exist for primes 
up to 11, but it is not clear that complete BES-MUBs depend 
in any way on the existence of such subgroups. The family of 
Cayley graphs under consideration (defined for all primes p) 
is actually the graph complement of a family of Ramanujan 
graphs, and we are able to list some general graph-theoretic 
properties that hold for all values of p. In sectionHJwe review 
the necessary background concerning the Clifford group and 
introduce some graph-theoretical concepts that will be useful 
in later sections. In section |lll| we explicitly give the recipe for 
constructing BES-MUBs and relate our work to a well-known 
MUB construction that uses finite field methods. Section iTVl 
further explores the quantities and concepts from graph theory 
that can be applied to our family of Cayley graphs, and finally, 
Appendix lAl provides a description of the MUB observables 
in terms of stabilizer measurements as well commuting sets of 
Pauli operators. 



The finite-dimensional analogues of position and momen- 
tum operators are denoted by X and Z, arbitrary products of 
which are called displacement operators D, indexed by a vec- 
tor it = (mi , U2) S Z 2 : 

X\j) = \j + 1) Z\j)=ri\j) (co = e 2 ™/") (1) 

D u = x" , " 2 X" l Z 1 ' 2 X = e (p+^IP. (2) 

The Weyl-Heisenberg group (or generalized Pauli group) for 
a single qupit is given by 



g p = {z c D u \ueZ 2 p ,ceZ p }. 



(3) 



The set of unitary operators that map the Pauli group onto it- 
self under conjugation is called the Clifford group (sometimes 
called the Jacobi group): 

c p = {ceu(p)\ug p tf = g p }. 

The fact that every Clifford operation in dimension p can be 
associated with a matrix F g SL(2,Z p ) in addition to a vector 
u £l? p results from the isomorphism 

Cp^SL{2,Z p )xZ% (4) 

established by Appleby ll32tl . where C is the Clifford group. If 
we specify the elements of F and u as 



I °! g J £SL{2,Z p ) u 



then Appleby provides an explicit description of the unitary 
matrix Ctp\ u \ E C p in terms of these elements i.e., 



C, 



(F\u) 



D U U F 



U F = { Vp 



(3 = 0. 



(6) 



(7) 



Note how composition and inverses can be represented in 
this notation 113111 



u+Fv) 

C {F\u) =C {F\u) =C (F- l \-F- l u) 



(8) 
(9) 



We will have need to relate the matrix trace Tr (C(p |«)) to 
the matrix trace Tr(F) modulo p: 



|Tr (Q 



{F\u)> 



e{0,^p,p} ifTr(F)=2 



(10) 



ll ifTr(F)^2 
To see why this is so we must define the Legendre Symbol 

{1 if x is a quadratic residue (mod p) 
— 1 if x is a quadratic non-residue (mod p) 
if x = (mod p). 
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and quote a result from Appleby ll32ll 

(Casel: p = 0^a^0) 

(\£ p (a)\ = l (Tr(F)^2) 
\Tt(C {f]u) )\= I \e p (y)\VPKo (Tr(F) = 2,y^0) 
[p5„ lj0 5 M2i o (Tr(F) = 2,y=0) 

(ID 

( Case 2: p ^ ) 

llWr \|JlW- 2)1 = 1 ( Tr (^)^ 2 ) 
1 lC(F| " )jl \|^(-P)lv^8» 2 ,p-i(i-«)» 1 (Tr(F)=2) 

(12) 

Finally, we note some important facts regarding the struc- 
ture of the group SL(2,Z p ). A minimal set of generators is 
e.g. 

SL(2,Z p ) = ^J J), (J j)^ (13) 

It has order \SL(2,Z p )\ = pip 2 - 1) and can be partitioned 
into p + 4 conjugacy classes 113311 . each of which has constant 
trace. If we partition SL(2,Z p ) by the matrix trace of its ele- 
ments, Tr(F), we see the following 



{F|£ p ((Tr(F)) 2 -4) 


= 1} 


p(p + l) 


(14) 


{F\i p ((Tr(F)) 2 -4) = 


-1} 


P(P-1) 


(15) 


{F\e p ((Tr(F)) 2 -4) 


= 0} 


p 2 


(16) 



The final sets {F|Tr(F) = 2} and {F\Tr(F) = -2} are each 
comprised of three conjugacy classes. Many of these facts 
will be used in subsequent sections, particulary section |IV] 
concerning graph-theoretical properties of Cayley graphs that 
are relevant to the construction of BES MUBs. 

B. Graphs: Cayley Graphs and Maximum Cliques 

We review some relevant notation and properties of graphs 
that can be found in any standard reference (e.g, lHH). An 
undirected Cayley graph Y(G,T) with an associated finite 
group G and set T C G, is the graph whose vertices are the el- 
ements of G and whose set of edges is {g\ ~ gi\gi g2 <= T}. 
We must have I ^ T and T = T. The resulting graph 
T(G,T) is regular i.e. each vertex has degree \T\, and the 
number of (undirected) edges is given by A|G||T| . A com- 
plete graph of order n, denoted K„, is a graph with n vertices, 
each of which is adjacent to every other vertex (see Fig[T](a) 
for an example K$). A subgraph, Y ', of Y, is a graph whose 
vertices form a subset of the vertices of Y and the adjacency 
relation is inherited from Y. A clique of Y is a complete sub- 
graph of r, where the size of the clique is given by the number 
of vertices in this subgraph. The largest possible clique (not 
necessarily unique) contained in Y is a maximum clique, the 
size of which is usually denoted C0(T). We discuss graph- 
theoretic properties, and what they say about the problem at 
hand, in more detail in Section. [IV] 



(IY, YI) 




(a) 



(b) 



FIG. 1. MUB structure: (a) MUBs in dimension 2 2 . Each box repre- 
sents a two-qubit stabilizer state p = i Y.seS s where S is the abelian 
subgroup generated by the Pauli operators contained in {•). Vary- 
ing the signs of the generators creates a complete orthonormal basis 
from each representative pair. Lines between boxes indicate that the 
overlap between two states is Tr(p a pi) = 4. The solid lines depict 
the complete graph on 5 vertices, A~5, and this corresponds to a com- 
plete MUB on this Hilbert space. The dashed lines depict a trian- 
gle, K 3 , which forms a BES MUB. (b) Two different BES MUBs 
(solid and dashed complete graphs K3) that partition 5L(2,Z2), 
where each 2x2 matrix F corresponds to the Jamiolkowski state 
(/(g)C (F | 0) )Ipo \jj)/V2. Adjacent vertices F U F 2 £ SL{2,Z p ) sat- 
isfy Ti(F~ l F2) 2, which in terms of the corresponding density 
matrices implies Tr^^p^) = ^. 



HI. CONSTRUCTION OF THE RESTRICTED MUB 

The goal is to create a set of states, S, of size |<S| = p 2 (p 2 — 
1 ) that is partitioned into p 2 — 1 subsets, where each subset, 
containing p 2 states, forms an orthonormal basis. Labeling 
the basis with a supercript and the individual states within a 
basis using a subscript we have 

5 = {| ¥ l)...|^)...| ¥ ;J- 1 )}. 
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This is a mutually unbiased basis if 



i<¥*i¥r 



1 



(1 -8t,n)+8jfe,„8; )/ 



The \yn) of the set S, that comprises our bipartite entangled 
stabilizer MUB (BES MUB), will be maximally entangled 
stabilizer states - formed by applying a Clifford operation, 
C, to one half of a maximally entangled state 

|/c) = (/®C) P £^. 

The overlap \(Jc„,\Jc„}\ between any two such states is 
given by 

l(/c m |/c n )| = i|Tr(4c„)|. 

Using the notation we have previously described, it is easy 
to show using Eq.s[8l— [TOlthat 



if Tr(F- l K) ^2 



then 



Tr 



C, 



{F\u)i 



c, 



(Jf|v) 



= 1 Vw,v£Z 2 



(17) 



i.e., a pair of matrices F,K £ SL(2,Z p ) satisfying 
Tt{F- 1 K)^2 defines a pair of mutually unbiased ba- 
sis. Since the subspace under consideration has dimension 
(p 2 - l) 2 , and since each basis contains p 2 — 1 independent 
states, we require a total of p 2 — 1 matrices Fj £ SL(2,Z p ), 
satisfying, pairwise, Tr(Ff l Fj) ^ 2, in order to create the 
BES MUB. 
Define 



G = SL(2,Z p ) 

T = {F€SL(2,Z p )\Tv(F)^2} 



\G\=p( P 2 -l) (18) 



then the Cayley graph T(G, T) has the property that two ver- 
tices Ft and Fj are adjacent if and only if Tr(F^~ l Fj) ^ 2. A 
clique of size p 2 — 1, if it exists, immediately gives the de- 
sired complete BES MUB by the preceding discussion. Fur- 
thermore, a clique of size p 2 — 1 must be a maximum clique 
since the dimension of the Hilbert space for local maximally 
mixed operators is (p 2 — l) 2 . 

Theorem 1 A pair of matrices F\,F2 £ SL(2,Z p ) satisfying 
Tr(F^ F2) 7^ 2 defines a pair of mutually unbiased bases in 
H 2 = C P (3C P (via the relationship between SL(2,Z p ) and 
the Clifford group). A set of matrices J- = {Ft}, of order 
\J-\ = p 2 — \, such that pairwise Tr(F^ l Fj) 7^ 2 (mod p), de- 
fines ( i) a complete bipartite entangled stabilizer MUB ( ii) a 
maximum clique of the Cayley graph defined in Eq. (fTS]). 



One can check using a computer algebra system 037L 
IHl] that the following subgroups H p < 5 , L(2,Z p ) have order 
\H P \ =p 2 -l, and every pair of elements Fi,Fj £ H p satisfies 



Tr(/^ Fj) 7^ 2 (i.e. these subgroups provide complete BES 
MUBs). 



p = 3: H 3 = 



p = 5: H 5 



p = 7: Hj 



p = ll: H n = 



oh f 1 1 



2 0/^12 
2 \ / 1 1 



2 J' \ 2 3 
2 \ / 1 1 



3 'U 5 



1 \ / 4 



10 / ' I 8 10 



(19) 
(20) 
(21) 
(22) 



In fact for every prime dimension p < 11 we can parti- 
tion SL(2,Z p ) by using p distinct max-cliques of size p 2 — 1. 
For odd primes it suffices to consider the left cosets of H p in 
SL(2,Z p ) where 



F> 



1 
t 1 



t gz b 



(23) 



are the left coset representatives. For p = 13 and higher, 
we were unable to find cliques saturating the upper bound of 
p 2 — 1. It is known that, for any primes p > 13, there does 
not exist a subgroup H p of size \H P \ = p 2 — 1 039I1 . but we are 
unaware of any proof that cliques of size p 2 — 1 (i.e. complete 
BES MUBs for p-dimensional systems) necessarily depend 
on this subgroup structure. A deterministic search for a clique 
of size 168 in the T(G, T) graph for p = 13 is infeasible, given 
the computational complexity of the max-clique problem. A 
heuristic search was able to find a clique of size 158, however. 

By adapting a well-known power-of -prime construction for 
complete MUBs (Bandyopadhyay et al. ITlll ) we can show 
that the size of the largest clique satisfies 



(0(T)>p( P -l) Vp 



(24) 



To be specific, Section 4.3.1 of Hi 1H describes the construc- 
tion of a complete set of MUBs for dimensions p 2 . In their 
notation, this amounts to finding a set of p 2 2 x 2 symmet- 
ric matrices {A} such that det(A ; - — Ak) 7^ 0. Suitable sets of 
matrices are parameterized by two elements s,t £ Z p via 



{A} 



{( 



-tb 



V a,b £ Z„ 



(25) 



A little thought reveals that every A with non-zero off di- 
agonal element b can be related one-to-one with a matrix 
F £ SL(2 : 1 p ), where F has a non-zero element (3 (F defined 
as per Eq. ©), 



F(a,b,s,t) 



-ab- 1 
b-a 2 b- l s- 



-1 



at —ab s — t 



One can check that the det(A ; — Ak) 7^ condition translates 
to Tr(Ff l Fi ( ) 7^ 2, as one would expect. In this way we can 
create a set of p(p — 1) matrices F £ SL(2,Z p ) that form a 
clique in our Cayley graph F(G,T). In general, sets of ma- 
trices formed this way cannot be extended with an additional 
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p — I matrices (having p^ = 0) to form a complete BES 
MUB i.e., they form (part of) a maximal, but not maximum, 
clique in T(G,T). However, we can often slightly improve 
upon the lower bound e.g., we can construct cliques of size 
p(p — 1) + 2 for primes up 17. A consequence of Eq. (l24l > is 
that the fraction of pairs (Fi,Fj) that do not define mutually 
unbiased bases, out of the total number of such pairs (F(,Fj), 
vanishes as p — > °°. In Appendix lAl we explicitly give the 
observables involved in these BES MUBs in terms of tensor 
products of Pauli operators. 



set is the largest such set (not necessarily unique) that can be 
found in the graph, and the independence number, oc(r), of 
a graph is the size of this maximum independent set. The 
discussion preceding Eq. (l26l i can equally well be interpreted 
as providing a lower bound on the independence number of 
r(G,r); there are p 2 — 1 independent sets of size p, wherein 
two elements F(, Fj that satisfy (a,-, p,) = (a,-, Pj) are pair- 
wise non-adjacent, hence 



a(r) > p Vp 



(27) 



IV. SOME GRAPH-THEORETIC PROPERTIES OF THESE 
CAYLEY GRAPHS 

In this section we further investigate the graph-theoretical 
properties of the family of Cayley graphs that were previously 
shown to be closely related to BES MUBs. Without loss of 
generality, the elements F 6 SL(2,Z p ) can be ordered lexico- 
graphically by the vectors constituting the rows of the matrix 
F i.e. 



{^} = {^i 



1 

-1 

-1 -1 
-1 -2 



1 

-1 1 



H("i:)}- 



It is easy to see that (i) there are p 2 — 1 possibilities for (a, P); 
(ii) each such (a, p) in turn allows for p possible (y, 8). Any 
two elements Fi, Fj, for which (a,-, P,) = (a,-, P ; ), cannot be 
connected by an edge since 



Tr 



a p 

Ji 5; 



-1 



a p 



det/v + detf) = 2. (26) 



The so-called vertex coloring problem for graphs involves as- 
signing a label (color) to every vertex of the graph, such that 
adjacent vertices cannot be assigned the same color. The min- 
imum number of colors required to do this is the chromatic 
number, denoted %(F). It is a basic fact l34ll that the chromatic 
number of a graph is bounded below by the clique number i.e. 
K>(r) < %(F). The discussion leading to Eq. [26l immediately 
implies that a p 2 — 1 coloring of the Cayley graph T(G, T) is 
possible: assign the same color to two vertices Fi, Fj if and 
only if (a,-, P/) = (a,-, P ; ). Since the chromatic number % is 
bounded below by the clique number C0(T), we know that this 
coloring is minimal for primes 2 to 11. Hence 



co(r) = x(r) = P 2 -i 
«(r)<x(r)<p 2 -i 



p€ {2,3,5,7,11} 
Vp 



Note that the upper bound co(r) < p 2 — 1 is a graph-theoretical 
inequality that confirms the geometrical argument preceding 
Theorem[T]i.e., the number of BES mutually unbiased bases 
that can fit in a Hilbert space Ti. p i = C p (g> C p is at most p 2 — 1 . 

A concept closely related to cliques and colorings is that of 
independence. An independent set of a graph is a set of ver- 
tices, no two of which are adjacent. A maximum independent 



The physical interpretation of this is that we can always find a 
set of p bases such that, pairwise, no two are mutually unbi- 
ased with respect to each other. 

The adjacency matrix of a graph Y with n vertices is an 
n x n matrix A[r] with elements Ay = 1 if vertices i and 
j are adjacent, and A,j = otherwise. Knowledge of the 
spectrum of an adjacency matrix often allows us to find, or 
bound, many quantities of interest. We denote the spectrum 
of the p(p 2 — 1) x p{p 2 — 1) adjacency matrices A[r(G,r)] 
as {A-g , A-'" 1 , A,™ 2 , A,™ 3 } where m,- denotes the multiplicity of 
A,-. The complement of a graph F, denoted F, is the graph 
with same vertex set as F, but where two vertices are adjacent 
in r if and only if they are not adjacent in F. The spectrum of 
a graph and its complement can be related in a simple way for 
the case of regular graphs (the case we deal with in this work), 
as the following theorem demonstrates. 



Theorem 2 (Brouwer and Haemers II36I1 ) Suppose F is a k- 
regular graph on n vertices with 4 distinct (adjacency) eigen- 
values {k = Ao > A; > A2 > A3}. If, in addition, both F and 
its complement, F, are connected, then F also has 4 distinct 
eigenvalues, {« — k — 1 > —A3 — 1 > — A2 — 1 > — Ai — 1}. 

The Cayley graphs we studied, defined in Eq. ( fT8l ). are ac- 
tually the graph complement of a well known family of graphs 
(that form a family of Ramanujan graphs, amongst other in- 
teresting properties), whose spectrum is known exactly. 

Theorem 3 (Lubotzky flU) Let G = SL(2,Z p ), and let T 
(i.e., the connection set for the Cayley graph) be the union 
of the conjugacy classes c\ and c v of the elements 



1 
1 1 



1 
V 1 



where V is a generator of the cyclic group Z* = Z p /{0}. Then 
T = {F e SL(2,Z p )\Tr(F) =2,F ^ I}, \T\ =p 2 -l and the 
spectrum of the corresponding Cayley graph A[r(G,7")] de- 
noted {A™ , A'" 1 , A™ 2 , A™ 3 } , is 



Ao — 

Ai = 

A 2 = 

A3 - 



P 2 -l, 

P-h 

0, 

-(P + 



m\ 

m 2 
;??3 



1 



= P 



(p-2)(p+l)72 

2 



p{p-\) 2 /2 
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Combining the two preceding theorems (connectedness is 
obviously satisfied by our Cayley graphs) allows us to com- 
pletely characterize the spectrum of the canonical Cayley 
graph Eq. (fT8l ) that we used to search for BES MUBs. 

Theorem 4 (Spectrum of graphs defined in Eq. ( 1781 )) Let 

G = SL(2,Z p ) and T = {F E SL(2,Z p )\Tr(F) ^ 2}. Then 
\T\ = \G\ — p 2 and the spectrum of A\T(G,T)] denoted 

h) = P(P 2 -1)-P 2 m = 1 

\\ = p mi = p(p-\) 2 /2 

A-2 = — 1 mi = p 2 

h = -p m 3 = {p-2)(p + l) 2 /2 

At this point we note that the problem of finding BES 
MUBs, framed as finding maximum cliques of size p 2 — 1 in 
the Cayley graph F defined by Eq. ( 1181 . is completely equiv- 
alent to finding maximum independent sets of size p 2 — 1 in 
the complement, T, of that graph i.e., 

3 complete BES MUB <^ co(T) = p 2 - 1 = a(T). 

Unfortunately, it seems that existing spectral lower bounds 
on the clique number are of little help for the task of prov- 
ing existence of BES MUBs. Nonetheless, using some well- 
known spectral bounds we list some implications for the 
graphs r(G, T) under consideration. A lower bound on the 
chromatic number is given by 

x(r)> 1-^=^-1), 

which, in conjunction with Eq. (126b . shows that 
p(p — 1) < %(r) < p 2 — 1. In fact, this lower bound 
was already implied by Eq. (|24| |. 

For a regular graph, F, on n vertices, Hoffman (unpub- 
lished) and Lovasz 14011 proved the formula 

. -nKun -nk-i 
a(r)<T — = r — =p + l, 

A-max — A- m i n Ao — A3 

which, in conjunction with Eq. d27b gives us 
p < Ot(r) < p + 1. As a final remark on spectral impli- 
cations, we note that the spectrum exhibited in Thm. [4] 
classifies F(G,T) as a so-called walk-regular graph fiUl . 

V. CONCLUSION 

We have shown how the set of bipartite entangled stabilizer 
(BES) states can be partitioned into sets of mutually unbiased 
bases (MUBs), whose span is sufficient and minimal to de- 
scribe an interesting class of operators that includes (mixtures 
of) Jamiolkowski states, Clifford witnesses 10] and more. Mu- 
tual unbiasedness of two stabilizer orthonormal bases is easily 
shown to be equivalent to a simple relation on pairs of ma- 
trices from SL(2,Z p ). Pairs of matrices satisfying this rela- 
tion are adjacent vertices on a naturally defined Cayley graph, 



and the problem of finding complete (optimal) BES MUBs 
is transformed into that of finding maximum cliques in the 
Cayley graph. In a different mathematical context, the graph 
complement of our Cayley graphs are well-studied, and so we 
can quote, for example, the exact spectrum of the adjacency 
matrix for all prime values p. The most interesting open ques- 
tion is whether such BES-MUBs exist for all primes, or in- 
deed for any primes greater than 1 1 . For the closely related 
task of finding minimal unitary designs, a discussion by Chau 
I2H1 (invoking Dickson's theorem on the existence of certain 
subgroups of SL(2, Z p )) suggests that minimal unitary designs 
only exist for primes up to 11. It remains to be seen whether 
the latitude afforded by seeking BES-MUBs, as opposed to 
subgroups of SL(2,Z p ), allows for construction of optimal 
BES-MUBs when p > 13. 

Appendix A: Measurement Operators for BES MUBs 

Given a matrix F E SL(2,Z p ), this defines an orthonormal 
basis J- in the bipartite Hilbert space H p i via 

T = {\j F u )yuez 2 p }, \(j f u \j*)\ = ?> u>v (ad 

where \J F U ) = (7® C (m ) £ M 

j=Q VP 

We will show how the basis T can be rewritten in terms of sta- 
bilizer measurements, and subsequently how T can be identi- 
fied as the simultaneous eigenbasis of a set of p 2 — 1 commut- 
ing Pauli operators. 

Using so-called symplectic notation, the general form for 
multi-particle stabilizer operators with vectors x = {x\ ,X2, ■ ■ ■ ) 
and z = (zi , Z2j • • • ) with x t , Zi E Z p is 

P(x\z) = (X Xl ®X X2 ...) (Z Z1 <g> Z Z2 . . . ) . (A2) 

Measuring a two-qupit Pauli operator corresponds to project- 
ing with a rank-p projector, II, 

n := ^-(x^ 2 \zi,z2)[k] = -{ I + (S) ' kp (x 1 ,x 2 \ zl ,z 2 ) +■■■ 

(A3) 

The product of two appropriately chosen such projectors, 
n,II', defines a rank-1 operator - a stabilizer state: 

l¥)(¥l = ^2 I> = nn '> 

p seS s 

where Q s = (g,g') is a subgroup, generated by two com- 
muting Pauli operators g and g\ of the group Q2 = 
{co 6 'P( A | 7 )|x,z E Z 2 ,c E Zp]. In symplectic notation g — 

( ti~ k P(x u x 2 \zi,z2) an d s' = 03 ~^ p (x! y \z' ,z' ) an d commutativity 
of g and g' reduces to 

^ XiZi - x\z\ = mod p. 

(=1,2 
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Given u — (111,112) € 1?„ and P 7^ 0, the following two sets 
of projectors are equal, up to re-ordering 

Vu- {\Ju)(Ju\} = {no.oiap-'.-p-^Miinfo.ii-p-'.p-'s)^]} 

When p = 0, the following two sets of projectors are equal, up 
to re-ordering 

Vm : {\J u ){ J u 1} = { n (l,a|0,7)[M 1 ] n (0,0[l,-8)[H 2 ]}- 

Many existing constructions for complete MUBs in "Hd 
(with power-of-prime dimension d) are based around the par- 
titioning of d 2 — 1 non-identity Pauli operators into d + 1 
classes, each of which contains d — 1 mutually commuting 
operators. Each basis within the MUB is then given by the si- 
multaneous eigenbasis of the d — 1 mutually commuting oper- 
ators (i.e., each class is associated with exactly one orthonor- 
mal basis, for a given partitioning). We can frame the con- 
struction of BES MUBs in this language too, with the modifi- 
cation that we are partitioning the set of all weight-two Pauli 



operators i.e. the subset {P {xi ^ ]z[ , Z2 )/{P[ XI , [ Z1 ,0) , P(o,x 2 \0,z 2 ) } } 
of size (p 2 — l) 2 . With individual classes containing p 2 — 1 
operators, there can only be at most p 2 — 1 such classes. It 
should be clear that a set of n matrices F £ SL(2,Z p ) (satisfy- 
ing Tr(F i ~ l Fj) 2) is equivalent to n non-overlapping classes 
of weight-two Pauli operators, each class containing p 2 — 1 
non-identity elements. Recalling Eq. ( IA1 b for the definition 
of the basis associated with F, then the associated class of 
unitary operators is the subgroup Q s = (g,g') of Q2- The si- 
multaneous eigenbasis of all p 2 Pauli operators in Q s forms an 
orthonormal basis. When P 7^ the class of Pauli operators 
corresponding to F is given by 

Gs(F) = (g,g') := (■P(i,0|ap-i,-p-')'' P (0,l|-p- I ,|3- 1 8))- 

When p = the class of Pauli operators corresponding to F is 
given by 

Gs(F) = (g,g') := (f(La|0j),-P(0,0|l,-8))- 
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